Well-developed deformation in 42 Si 
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Excited states in 38 > 40 > 42 Si nuclei have been studied via in-beam 7-ray spectroscopy with multi- 
nucleon removal reactions. Intense radioactive beams of 40 S and 44 S provided at the new facility of 
the RIKEN Radioactive Isotope Beam Factory enabled 7-7 coincidence measurements. A prominent 
7 line observed with an energy of 742(8) keV in 42 Si confirms the 2 + state reported in an earlier 
study. Among the 7 lines observed in coincidence with the 2 + — > + transition, the most probable 
candidate for the transition from the yrast 4 + state was identified, leading to a 4^ energy of 2173(14) 
keV. The energy ratio of 2.93(5) between the 2^~ and 4^ states indicates well-developed deformation 
in 42 Si at N = 28 and Z = 14. Also for 38 - 40 Si energy ratios with values of 2.09(5) and 2.56(5) were 



obtained. Together with the ratio for 42 Si, 
isotopes from N = 24 to N = 28. 



the results show a rapid deformation development of Si 



PACS numbers: 25.60.-t, 23.20.Lv, 27.40.+Z, 29.38.Db 



Shell closures and collectivity are important proper- 
ties that characterize the atomic nucleus. Interchange 
of their dominance along isotopic or isotonic chains has 
attracted much attention. The recent extension of the 
research frontier to nuclei far away from the valley of sta- 
bility has revealed several new phenomena for neutron- 
or proton-number dependent nuclear structure. For ex- 
ample, a weakening or even disappearance of shell clo- 
sures occur in several neutron-rich nuclei at N = 8 [3-0| 
and N = 20 A well known example in the case of 

N = 20 is the so-called 'island of inversion' Q located 
around the neutron- rich nucleus 32 Mg. The low excita- 
tion energy of the first 2 + state E x {2^) and large E2 
transition probability clearly indicate shell quench- 
ing in 32 Mg despite the fact that N = 20 is traditionally a 
magic number. The next magic number, TV = 28, which 
appears due to the ,£7/2-/5/2 spin-orbit splitting, has also 
been explored [8l4l3|. Weakening of the shell closure is 
seen by the decrease of the 2+ energy for N = 28 iso- 



tones from 3.83 MeV in the doubly-closed nucleus 48 Ca 
to 1.33 MeV in 44 S as the proton number decreases from 
Z = 20 to Z = 16. Another region of shell stability 
has been shown to exist at N, Z = 14 due to the d 5 / 2 - 
d 3 / 2 spin-orbit splitting. The silicon isotopes from 28 Si to 
34 Si exhibit relatively high l\ energies ranging from 1.78 
MeV ( 28 Si) to 3.33 MeV ( 34 Si) reflecting the Z = 14 sub- 
shell closure. However, the 2+ energy gradually decreases 
from 36 Si to 40 Si [lj|, suggesting a development of 
quadrupole collectivity for isotopes with N > 20. 

With proton number Z = 14 and neutron number 
N = 28, the nuclear structure of 42 Si is of special interest. 
A simple but important question that arises is whether 
the weakening of the N — 28 shell closure continues, 
causing an enhancement of nuclear collectivity, or if shell 
stability is restored owing to a possible doubly magic 
structure. A study on Si was made by a two-proton re- 
moval reaction experiment with radioactive 44 S beams at 
the NSCL [lij]. The small two-proton removal cross sec- 
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tion was interpreted as evidence for a large Z = 14 sub- 
shell gap at N = 28, indicating a nearly spherical shape 
and a doubly closed-shell structure for 42 Si. Contrary to 
this result, a disappearance of the N — 28 spherical shell 
closure around 42 Si was concluded by another experimen- 
tal study at GANIL with the same reaction [l|| owing to 
the observation of a low-energy 7 line interpreted as the 
transition from the 2+ state at 770(19) keV to the ground 
+ state (0+ s ). This low E x (2^) value supports the non- 
magic nature of 42 Si expected from comparison of the 
/3-decay half-life of 42 Si with QRPA calculations [ll|. 

The disappearance of the N = 28 shell closure for 
42 Si was theoretically pointed out in several recent stud- 
ies with shell-model and mean-field approaches 17- 19$ . 
These studies predict a well-developed large deformation 
of the ground and low-lying states. For further under- 
standing of the structure of 42 Si and, more generally, the 
mechanism of interchange between the shell closures and 
quadrupole collectivity along the N = 28 and Z = 14 
chains, further experimental input is necessary. In ad- 
dition to the energy of the 2+ state, the location of the 
4^~ state is of crucial importance, since the energy ratio 
between the 4^" and 2^ states (-R4/2) represents the char- 
acter of quadrupole collectivity: a ratio of 2 is expected 
for harmonic vibration and 3.33 for rigid-body rotation, 
as extremes. Hence, the systematic study of E x (2j) and 
E x (^i) values is useful to deepen our understanding on 
the mechanism for the evolution of nuclear structure in 
the vicinity of 42 Si. 

The purpose of the present study is to find the 4j~ state 



in Si as well as to revisit the 2{ 



transition. Ow- 



ing to the secondary 44 S beam with the world-highest 
intensity provided by the RIKEN Radioactive Isotope 
Beam Factory (RIBF), population of the 4+ state in 42 Si 
by the two-proton removal reaction was possible, and 
even a 7-7 coincidence analysis was enabled to establish 
the level scheme with the help of the high- efficiency detec- 
tor array DALI2 (Detector Array for Low-Intensity radi- 
ation 2) 2(1 2l|. Additionally, we studied multi-nucleon 
removal reactions of 44 S and 40 S to populate 4 + states in 
40 Si and 38 Si in order to obtain valuable information on 
the evolution of the quadrupole collectivity in neutron- 
rich Si isotopes provided by the systematic trend of the 
ratio i?4/2- 

The experiment was performed at the RIBF operated 
by RIKEN Nishina Center and the Center for Nuclear 
Study, University of Tokyo. A primary 48 Ca beam at 
345 MeV/nucleon with an average intensity of 70 pnA 
bombarded a 15-mm-thick rotating beryllium target. A 
secondary 44 S or 40 S beam was produced by projectile 
fragmentation and analyzed by the BigRIPS fragment 
separator [22| as in earlier experiments [2^, [24[. The 
energy and intensity of the secondary 44 S ( 40 S) beam 
was approximately 210 MeV/nucleon (210 MeV/nucleon) 
and 4 x 10 4 particles per second (pps) (6 x 10 4 pps), 
respectively. 
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FIG. 1. (a) Doppler-shift corrected 7-ray spectrum obtained 
for the C( 44 S, 42 Si 7) reaction, (b) 7-7 coincidence spectrum 
gated with 742 keV 7 ray. Solid curves in both figures show 
the results of fits using Gaussian functions (dashed curves) 
and an exponential curve as the background (dotted curve) by 
the method % 2 minimization. The dot-dashed line indicates 
the expected line shape in the case that the 2032-keV 7 ray 
decays solely to the 2 + state. 



A carbon target with a thickness of 2.54 g/cm 2 was 
located at the F8 focus for secondary reactions. Reac- 
tion products were analyzed using the ZeroDcgree spec- 
trometer 



25], and identified using the energy loss (AE), 



magnetic rigidity {Bp), and time-of-fhght (TOF) infor- 
mation. The Bp value was obtained from the position 
at the dispersive focus F9 measured by parallel plate 
avalanche counters (PPACs) [26]. AE was measured by 
an ionization chamber at the achromatic focus Fll and 
the TOF information was obtained from the time differ- 
ence between plastic scintillators at F8 and Fll. The 
inclusive cross section for the C( 44 S, 42 Si) reaction at 210 
MeV/nucleon was obtained to be 0.15(2) mb, which in- 
dicates a monotonic rise of the cross section as compared 
with 0.12(2) mb at 98.6 MeV/nucleon and 0.08(1) 
mb at 39 MeV/nucleon [l6j]. 

De-excitation 7 rays were detected by the DALE ar- 
ray in coincidence with the outgoing 42 Si, 40 Si, and 38 Si 
particles. DALE consists of 186 Nal(Tl) detectors sur- 
rounding the reaction target that cover an angular range 
of 11° — 160° with respect to the beam axis. Typical 
photo-peak efficiency and energy resolution are 20% and 
10% (FWHM), respectively, for 1 MeV 7 rays emitted 
from nuclei moving with /3(= v/c) ~ 0.6. These val- 
ues were obtained by Monte-Carlo simulations using the 
GEANT4 code [H. 
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Figure [Ha) shows the 7-ray spectrum for the 
C( 44 S, 42 Si 7) reaction, where Doppler-shift effects have 
been corrected for. For each 7-ray peak, the energy was 
obtained by a fit of a Gaussian function and an expo- 
nential background curve, where the energy resolution 
was fixed to the simulated value. As seen in the figure, 
a predominant peak is observed at 742(8) keV, where 
the error includes statistical and systematic uncertain- 
ties. The systematic error is attributed to the uncertain- 
ties in the energy calibration (3 keV), the Doppler-shift 
correction (3 keV), and the ambiguity caused by a delay 
of 7-emission resulting in a shift of the source position 
(estimated to be 2 keV as a maximum lifetime of ~ 40 ps 
based on systematics [IH])- The observed peak energy 
agrees within 1.5 standard deviation with the value of 
770(19) keV reported in the study at GANIL pj|. In the 
higher energy region, three additional 7 lines have been 
identified in the present study at 1431(11), 2032(9), and 
2357(15) keV. 

In order to identify the transitions feeding the 2+ state, 
a 7-7 coincidence analysis was performed. The spectrum 
obtained from a gate on the 742-keV line is shown in 
Fig. [IJb). A clear peak is seen at 1431 keV, which may 
feed the 2f state directly from a higher-lying excited 
state at 2173(14) keV. By considering the 7-ray detec- 
tion efficiency, the yield of the peak is consistent with 
a 100% feeding of the 2+ state. A peak-like structure 
around E 1 ~ 2 MeV could correspond to the 2032-keV 
7 line observed in Fig. [lja) , providing another candidate 
that populates the 2 4 state. However, the yield of the 
peak in the 7-7 spectrum is lower than the expected value 
(indicated by the dot-dashed line), based on the inten- 
sity measured in Fig. QTa) and assuming a decay branch 
of 100% to the level at 742 keV. This suggests that the 
2032-keV 7 ray does not, or at least does not fully, pop- 
ulate the 2+ state. From the known tendency that yrast 
states, including the 2^~ and 4^~ states, are preferentially 
observed with larger cross sections in nucleon removal re- 
actions 29p 31 1 , those two 7 lines are possible candidates 
for the 4^ — > 2^ transition in 42 Si. Since the 7-7 coinci- 
dence analysis indicates that the 1431-keV 7 ray directly 
feeds the 2+ state as discussed above, 2173 (742 + 1431) 
keV for the 4+ energy is more probable among the two 
possibilities. Thus, we tentatively assign the state at 
2173(14) keV. The resultant R 4/2 value of 2.93(5) for 42 Si 
is rather close to the rigid-rotor limit. This contradicts 
the possibility of a doubly-closed structure suggested by 
the two magic numbers Z — 14 and N — 28, but sup- 
ports enhanced quadrupole collectivity in Si expected 
from the measured 2^~ energy [l|| and theoretical calcula- 
tions p~7i - figj . Furthermore, the large R 4 / 2 value indicates 
a large static ground state deformation of 42 Si. 

A search for 4^ states in 38 > 40 Si was also conducted. 
Figures [2ja) and (b) show the Doppler-corrected spec- 
tra for 38 Si and 40 Si, measured with the C( 40 S, 38 Si 7) 
and C( 44 S, 40 Si 7) reactions, respectively. In the case of 
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FIG. 2. Doppler-shift corrected 7-ray spectrum obtained in 
(a) the C( 40 S, 38 Si 7) reaction and (b) the C( 44 S, 40 Si 7) reac- 
tion. 



38 Si, the major peak likely consists of three unresolved 
7-ray lines, according to the ones observed in a previ- 
ous experiment 0. Their energies, 1071(12), 1168(22), 
and 1284(26) keV, were obtained by fitting the spectrum 
shown in Fig. [Ha), where three peaks with fixed widths 
and initial centroid positions estimated from Ref. [l3[ 
were used in the fitting procedure. The 1071-keV 7 line 
corresponds to the known 2^~ — > 0^~ s transition, while the 
1168- and 1284-keV lines are candidates for the 4^ — > 2^ 
7 ray. Here, we assign the 2239-keV state to be the most 
probable candidate for the 4+ state, based on the same 
arguments of yrast feeding that were discussed for 42 Si. 
We note that the alternative if assignment to the 2355- 
keV state does not affect the discussion on the systemat- 
ics of the level energies of Si isotopes discussed later. 

For 40 Si, three 7 lines are seen at 629(8), 985(11), and 
1539(16) keV in Fig. C^b). The 629- and 985-keV lines 
have been observed in the p( 42 P, 40 Si+7) reaction. The 



985-keV 7 ray was assigned to the 2{ — > Qj s transi- 
tion [lj|. The line at 1539 keV is reported here for the 
first time. According to 7-7 coincidence analysis, the 
629- and 1539-keV 7 lines are considered to be transi- 
tions from the excited states at 1614(14) and 2524(19) 
keV to the l\ state. Using similar arguments to those 
given for 38 Si and 42 Si regarding the preferential popula- 
tion of yrast states, either of the two states could be the 
yrast 4 + level. However, systematic trends of 2 + and 4 + 
energies suggest that the level at 2524 keV is the more 
likely of the two. The 1614-keV state is consistent with 
the previously observed tentative 1624(7)-keV state [l2j. 
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FIG. 3. (a) Ratio between the energies of the l\ and 4+ 
states (-R4/2) for Si isotopes. The horizontal lines at 2.0 and 
3.3 indicate the vibrational and rotational limits, respectively, 
(b) Excitation energies for 2+ and 4^ states, which are indi- 
cated by circles and squares, respectively. Filled symbols are 
results of the present study, and solid and dashed lines repre- 
sent predictions of the SM with SDPF-MU [H and SM with 
SDPF-U-MIX [H, respectively (see text for details). The 2+ 
energies of the N = 24, 26. 28 Si isotopes have been measured 
in previous works [2, 13 13 ■ 



On the basis of the large-scale shell model (SM ) calcula- 
tions in a Tr(sd) z ~ 8 v(pf) N ~ 20 model space [12j this level 
could be either a + or 2 + level. 

Figure [3] shows the isotopic dependence of the exci- 
tation energies of the 2^ and 4^ states together with 
their ratio R4/2 for 36-42 Si, where filled symbols repre- 
sent the present results, and open symbols for 36 Si are 
taken from Ref. [34j |. As seen in the figure, the ratios for 
36 Si and 38 Si [R 4/2 = 2.024(4) and 2.09(5)] are close to 
the vibrational limit, suggesting a nearly spherical shape, 
whereas R4/2 for 



40 Si increases to 2.56(5), indicating a 
deviation from the spherical shape or enhancement of 
quadrupole collectivity at N — 26. The lowering of the 
2^ energy as well as the increase of the reduced transi- 
tion probability or the deformation parameter obtained 
for 36 > 38 > 40 Si was interpreted as a narrowing of the N = 28 
shell gap [l^-Q] • The prediction of the SM calculation 
using the recent SDPF-U-MIX interaction [32j], an up- 
dated version of SDPF-U [l8| , is indicated by the dashed 
lines in Fig. [3j The new interaction, which allows np- 
nh excitations across the N = 20 shell gap, reproduces 
4/2 in a satisfactory manner up to 40 Si, but 
then deviates significantly from the experimental result 
for 42 Si in Fig. GHa). In the case of the N — 28 isotope, 
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FIG. 4. (Color online) Summary of the experimental and the- 
oretical levels in 38,40,42 Si 

isotopes. The dotted lines with the 
label S n indicate the location of the experimental (for 38,40 Si) 
or evaluated (for 42 Si) neutron separation energy [33]. Among 
the experimental levels, the lines with red color indicate the 
new results from the present study. The labels A, B, C, and 
D represents res pec tively the results of the shell model calcu- 
lations of Refs. [13, [H and [H, and mean field calculation 
of Ref. GJl. 



E x (2f) is lower than 40 Si and R4/2 further increases to 
2.93(5) despite the neutron magic number TV = 28. This 
indicates that a development of nuclear deformation con- 
tinues up to at least N — 28. In addition, the quadrupole 
collectivity increase in proton deficient N — 28 isotones 
turns out to continue to 42 Si with Z — 14. Thus, for 
42 Si well developed deformation has been experimentally 
established, and the possibility of a doubly-magic nature 
has been excluded. 

The solid lines in Fig.[3]indicate the results obtained by 
the SM calculation using the SDPF-MU interaction (l7j . 
which includes the tensor force in the effective interac- 
tion. The model reproduces the overall trends well, par- 
ticularly for the 4+ energies, where better agreement is 
achieved compared to SDPF-U-MIX. In particular, it pre- 
dicts a large R4/2 value of 2.8 for 42 Si which is close to 
the experimental result of 2.93(5). The mean field cal- 
culation with the relativistic energy density functional 
DD-PC1 [19( predicts a rotational band in 42 Si. Though 
the excitation energies of the 2^ and 4^~ states are about 
two times larger than the experimental values, the R4/2 
ratio (2.7) is not far from the present result. 

Figure |4] summarizes the experimental and theoretical 
level-energies for the isotopes 38 ' 40 ' 42 Si. The dashed lines 



indicate the levels observed in earlier studies [12], [l3j and 
the red lines show the 4 + states tentatively assigned in 
the present study. The labels A, B, C, and D represent 
theoretical calculations by the SM with SDPF-MU fl7jj . 



SM in Ref. [l2j . and mean 



SM with SDPF-U-MIX [: 

field calculation [lij], respectively. As seen in the figures, 
some model calculations predict more levels compared 
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with the ones experimentally observed. The experimen- 
tal energies of the calculated states with no empirical 
counterparts may provide a deeper understanding of nu- 
clear structure. Together with possible measurements on 
7-ray angular distributions and/or correlations, further 
efforts to identify those states present one target for fu- 
ture experiments with higher statistics by improvement 
of RI beam production. Other challenges on the exper- 
imental front, such as measuring excited states in 40 Mg 
and 44 Si, are encouraged to further trace the quadrupole- 
collectivity development along the chains of N — 28 iso- 
tones and Z = 14 isotopes. 

In summary, excited states in 38,40,42 nave been in- 
vestigated via in-beam 7-ray spectroscopy with multi- 
nucleon removal reactions in inverse kinematics by using 
210-MeV/nucleon 40 < 44 S beams. With the high-efficiency 
detector array DALI2 and the high intensity secondary 
beams provided at RIKEN RIBF, measurements with 
high statistics were achieved. The energy of the first 
2 + state in 42 Si was measured to be 742(8) keV, and the 
most probable candidate of the 4^ state was found at 
2173(14) keV with the aid of a 7-7 coincidence analysis. 
The 4+ states in 38 Si and 40 Si were assigned excitation 
energies of 2239(25) and 2524(19) keV, respectively. The 
systematics of the ratio i?4/ 2 of the 4+ and 2^~ energies 
in silicon isotopes from N — 24 to N — 28 shows a rapid 
development of deformation. The R4/2 value of 2.93(5) 
for 42 Si is the largest also among the known N = 28 
isotones, indicating that this nucleus has a character of 
well-deformed rotor despite the magic numbers N — 28 
and Z = 14. 
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